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Abstract. Radio observations at multiple frequencies have detected a significant isotropic 
emission component between 22 MHz and 10 GHz, commonly termed the ARCADE-2 Excess. 
The origin of this radio emission is unknown, as the intensity, spectrum and isotropy of 
the signal are difficult to model with either traditional astrophysical mechanisms or novel 
physics such as dark matter annihilation. We posit a new model capable of explaining the 
key components of the excess radio emission. Specifically, we show that the re-acceleration of 
non-thermal electrons via turbulence in merging galaxy clusters are capable of explaining the 
intensity, spectrum, and isotropy of the ARCADE-2 data. We examine the parameter spaces 
of cluster re-acceleration, magnetic field, and merger rate, finding that the radio excess can 
be reproduced assuming reasonable assumptions for each. Finally, we point out that future 
observations will definitively confirm or rule-out the contribution of cluster mergers to the 
isotropic radio background. 
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1 Introduction 

Observations of the extragalactic radio background have detected an isotropic excess in the 

frequency range of 22 MHz to 10 GHz [1-7]. The aggregate data, assembled by the Absolute 

Radiometer for Cosmology, Astrophysics, and Diffuse Emission-2 (ARCADE-2) Collabora¬ 

tion, show that the emission has a spectrum that falls as T~ 1 2 * * * ' 6 [6]. The overall intensity of 
the emission (hereafter referred to as the ARCADE-2 excess), exceeds the intensity expected 
from the isotropic portion of the galactic diffuse emission by a factor of a few [5, 8]. 

Several models have been posited in order to explain the intensity of the ARCADE-2 signal, 
including contributions from radio galaxies, radio supernovae, radio-quiet quasars, and dif¬ 
fuse electrons in the inter-galactic medium [9]. However, the extrapolation of these source 
classes to low-luminosities has failed to reproduce the total intensity of the ARCADE-2 ex¬ 
cess while maintaining consistency with constraints from infrared and X-Ray emission [10]. 
In addition to astrophysical sources, dark matter models have also been formulated in order 
to explain the ARCADE-2 excess. Dark matter models are well motivated for explanations 
of ARCADE-2 since they naturally produce large fluxes of relativistic e + e~ while producing 
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negligible thermal emission. This allows dark matter induced emission to simultaneously fit 
the intensity of the ARCADE-2 excess while remaining consistent with infrared constraints, 
and several similar models have been formulated to date [11-14], 

However, a second peculiar feature of the ARCADE-2 excess pertains to the high isotropy 
of the extragalactic signal. It was first noted by Holder [15] that the high isotropy of the 
ARCADE-2 excess may be in tension with the small total-anisotropy observed in high angular 
resolution radio observations of the extragalactic sky. In fact, even making the conservative 
assumption that the ARCADE-2 emission dominates the total anisotropy observed in VLA 
and ACTA observations, Holder placed an upper limit on the ARCADE-2 anisotropy which 
falls below level expected for emission sources that trace large scale structure [15]. Addition¬ 
ally, an analysis by Vernstrom et al. [16] calculated the luminosity distribution of detected 
radio sources at an frequency of 1.75 GHz, and found that a population of point sources 
would be unable to explain the total intensity of the ARCADE-2 emission unless they have 
individual emission intensities smaller than 1 /rJy. These observations strongly rule out most 
astrophysical, and dark matter source classes, requiring an extremely large number of very 
faint emission sources in order to produce the ARCADE-2 signal. 

One caveat common to both of the anisotropy constraints concerns the small angular held 
of view of radio interferometers. Due to the zero point subtraction necessary for this in¬ 
strumentation, they are insensitive to emission which is isotropic on larger angular scales. 
Specifically, [15] constraints the total anisotropy level on angular scales ranging from l = 6000 
- 28000, while [16] quotes a sensitivity to objects smaller than 2 arcminutes. Notably, this 
angular scale is smaller than that of the most massive galaxy clusters. In the context of 
dark matter models for the ARCADE-2 excess, it was noted by [14] that the anisotropy 
from this emission may fall below the isotropy constraints of [15] and [16] in scenarios where 
synchrotron radiation from dark matter annihilation is dominated by galaxy clusters, with 
significant emission extending to several times the virial radius of these massive objects. 
While this model is intriguing, it requires large magnetic held extensions in cluster sized 
objects, compared to typical cluster models [17]. In addition, this model also requires large 
dark matter substructure boost factors, which is consistent with the results of the Aquarius 
simulation [18-20], but is disfavored by [21]. 

An alternative model, first noted in [14], concerns the re-acceleration of non-thermal elec¬ 
trons through hydrodynamic shocks produced in cluster mergers. In the context of the dark 
matter models originally considered by [14], these shocks were employed to re-accelerate 
electrons originally produced via dark matter annihilation, significantly amplifying the total 
synchrotron flux produced via a dark matter annihilation event. This mechanism is capa¬ 
ble of smearing out the resulting synchrotron signal via re-acceleration, decreasing the total 
anisotropy of the synchrotron emission. In this paper, we show that cluster mergers can 
independently produce the intensity, spectrum and morphology of the ARCADE-2 excess 
without the need for seed-electrons from dark matter annihilations. In this case, we note 
that any population of non-thermal seed electrons can be efficiently re-accelerated to GeV 
energies. Moreover, the cluster mergers also produce significant magnetic fields within the 
cluster, allowing the re-born electrons to produce significant synchrotron radiation in the 
ARCADE-2 band. Finally, since strong mergers are most prevalent in the most massive 
clusters, the anisotropy from this emission mechanism occurs primarily on scales larger than 
the anisotropy studies of [15, 16]. Indeed turbulent reacceleration is one of the most popular 
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models employed to explain the observed radio emission at cluster scales [22], 

The outline of the paper is as follows. In Section 2 we discuss the physical model that produces 
synchrotron emission in cluster mergers, breaking down the process into (i) production of 
hydrodynamic turbulence and Alfven waves, (ii) the Alfven re-acceleration of e _ and e + in 
the intra-cluster medium, (iii) profiles of the density and the magnetic fields of clusters and 
the (iv) cluster merger rate. In Section 3 we compare our resulting model for synchrotron 
radiation against the observable parameters of the ARCADE-2 excess. Finally in Section 4, 
we discuss several caveats to our analysis and conclude. 


2 Models 

The production of the synchrotron emission in cluster mergers proceeds through several steps. 
Cluster mergers first induce fluid turbulence which then produces Alfven waves. These Alfven 
waves accelerate existing non-thermal electrons to GeV energies. The fluid turbulence also 
amplifies existing cluster magnetic fields, causing the GeV electrons to efficiently lose energy 
to synchrotron radiation in the GHz band. We will discuss the assumptions and outputs of 
our model in each step of the synchrotron production process. 

2.1 Turbulence and Alfvenic waves 

In order to calculate the production of Alfven Waves via fluid turbulence in cluster shocks, we 
closely follow the resonant acceleration model studied in [23] and [24]. Specifically, we assume 
that the energy spectrum of the turbulence induced by a cluster merger can be described by 
a power law: 

( \ —m 

KoJ (2 ' 1} 

where Wf is the energy density of the turbulent fluid at a wavenumber re, where re = 2ir/l is 
the wavenumber for a length scale l. The value kq = 2tt/Iq corresponds to the largest eddy 
size of the system, Iq. The energy density of the turbulence is given by Et ~ pv ^, where p is 
the mass density of the fluid and vt is the turbulent velocity at the largest scale Iq. Noting 
that the total energy density is the sum of the scales re, we have E t = jWf{n)dK and thus: 

W° f = (m - 1 )pvt/K 0 (2.2) 

Note that in this calculation we have ignored the upper limit for re as most of turbulent energy 
resides in the largest scale. The index of the turbulence, m is set to 5/3 for Kolmogorov 
turbulence and 3/2 for Kraichnan turbulence [25]. 

This fluid turbulence generates Alfven waves through Lighthill radiation [26, 27] . Continuing 
to follow [23] and [24], we limit the scale for the largest Alfven waves to correspond to the 
Taylor length, which defines the transition between large-scale ordered fluid motions and 
small-scale disordered motions: 


It 



(2.3) 
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Correspondingly, kt = 2 tt/It defines the minimum wavenumber of the Alfven mode. Here 
Re = lovt/v*. is the Reynolds number of the ICM, where u K = u p A e ff/3 is the kinetic vis¬ 
cosity of ICM; rripUp /2 = ksT is the thermal velocity of protons, where T is the ICM 
temperature. The value A e g = A 2 /A c is the mean free path for the transverse drift of 
protons in cluster magnetic fields, where X g = u p m p c/eB is the proton gyroradius and 
A c = 3 3 / 2 (/cbT) 2 /(47t 1 / 2 n p e 4 In A) is the mean free path for Coulomb collisions [28], where 
In A « ln(127rn p A£) e ) is the Coulomb logarithm and ADe = ^/ksT/n p e 2 is the Debye length. 
For typical values of the above parameters we obtain a Reynolds number: 

Re - iom (^y (soo^ft) m < 2 - 4 > 

x(icF^r(^v ) ,/2 (^) 2 


For each wavenumber k exceeding the threshold kt, Alfven waves are generated at a wavenum¬ 
ber k, calculated as in [24]: 


k = k 


v f (n) 


va 


(2.5) 


where Vf(n) is the velocity of an Eddy with wavenumber n, which can be calculated by 
Vf ~ {kW f(n)/p) 1 ^ 2 . Plugging in Equations 2.1 and 2.2 we get 


«/(«) = vt(m- 1) 1/2 


(l—m)/2 


and va is the Alfven velocity in plasma with magnetic field B [29], 
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Assuming that the energy injected into Alfven waves from turbulence follows a power law 
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where I a has the units of energy per volume per time per wave number. Within the context 
of Lighthill theory, St can be calculated as 


st 


3 [m — 1) 
3 — m 


which, combined with Equations 2.1 and 2.2, gives [23, 24]: 


h = r)A(s t - 1 )pv\ 



3 /( 3 —m) 
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where R = (tn — l)^ 1 (hq/kt) 1 ^ 171 , and ija ~ |2(3m — 4)/(3 — m)\ is a factor of order unity [30]. 
The portion of the total turbulent power going into the Alfven waves can be calculated by 
Pa = Jt:r lA<ik = pa ( v t / v 'aR) Pv\kt■ Assuming the same parameters as in Equation 2.4, 
the power is calculated as: 
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P A = 4.2 x 10 
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As a comparison, assuming that the power injection takes ~ 1 Gyr, the Alfven energy density 
is just 0.04% of the thermal energy density, £ t h = 3.2 x 10 -12 (n p / 10~ 3 cm _3 )(T/2 keV) erg cm -3 . 

The spectrum of the Alfven waves, Wk(t), responds to this energy input, and evolves due to 
both wave-wave and wave-particle coupling. These processes can be described by a continuity 
equation [31]: 


= ~j2 nw k (t)+i A (k , t) (2.i2) 

2—1 

The first term denotes the wave damping with relativistic particles in the ICM, which we 
will discuss next in Sec. 2.2. We note that in a complete description of the system, a cascade 
term d [DkkdWk(t)/dk\ /dk should also be included on the right-hand side of Equation 2.12. 
This physically corresponds to a wave-wave interaction term, with Dkk being the wave-wave 
diffusion coefficient. In the remainder of this work, we follow [24] and ignore this term in our 
analytic approximation. However, this approximation depends on both the energy density 
of the waves and the relativistic particle spectrum. We refer readers to [30] for a detailed 
comparison of the relevant time-scales over a number of physical environments. 

2.2 Alfvenic Acceleration 

Alfven waves can accelerate particles via resonant damping. In order to resonate with Alfven 
waves, an input population of non-thermal electrons is needed [32]. While we do not com¬ 
ment specifically on the origin of this charged particle population, we note that non-thermal 
protons and electrons are known to populate galaxy clusters. Possible astrophysical injec¬ 
tion mechanisms for these relativistic particles include collisionless shocks in merger events 
[33, 34], active galactic nuclei [35] and powerful galactic winds during starburst activity [36] 
(We refer readers to [22] for a detailed review). In our model, we simply assume that there is 
a sufficient density of non-thermal electrons to be re-accelerated by the turbulence powered 
by major mergers. 

The resonance condition between a wave with wave number k and a particle with mass m 
and momentum p can be written as [31] 


k 


flm 1 
P ppv A /c 
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where fl = eB/me is the cyclotron frequency, p = v\\/c is the projected particle velocity along 
the magnetic field, and the minus (plus) sign corresponds to electrons (protons) as damping 
particles. Notice that we have used the dispersion relation oj ~ k\\VA and assumed k « fen, 
which is reasonable for Alfven waves in the ICM [30]. 

The resonance condition indicates a relation between the minimum wave number and the 
maximum Lorentz factor of electrons. Considering that va c and \p\ is of order unity, 
Equation 2.13 implies 7 max ~ Op/ Lin = fl e /kT- From the dehnitions of Equations 2.3 
and 2.5 we obtain: 


1 nv A /Re\ (m_3)/4 
^ max (m — l ) 1 / 2 kq evt V 15 / 


(2.14) 


assuming consistent physical parameters with Section 2.1, this corresponds to a maximum 
electron energy: 
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(2.15) 


i? max represents the highest energy of an electron that could resonate with the turbulence, 
but can be considered as an estimation of the maximum energy of accelerated electrons as 
we will demonstrate at the end of this section. 


Using the resonance condition, the damping rate T(k) in Equation 2.12 due to the interaction 
with electrons can be written as [24, 37]: 


m = 


4ir 3 e 2 v\ 

c 2 k 
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(2.16) 


Following the argument in Section 2.1, the wave-wave interaction term in Equation 2.12 can 
be ignored, so the electron number density and the Alfven wave spectrum evolve according 
to the following coupled equations [23]: 


dW k (t) 

dt 

d_l 

dt 


-r (k)w k (t) + i A (k,t) 


ld_ 

p dp 



dj_ 

dp 


+ Sp 4 f 


(2.17) 

(2.18) 


where f(p) = N(p)/(Airp 2 ) is the phase space electron distribution; Sp 2 c corresponds to 
the synchrotron and inverse Compton emission power of an electron with S = 4 (B 2 + 
^CMB) e4 /(9mfc 6 ); -Bcmb = 3.25(1 + z) 2 p G is the equivalent magnetic field strength of 
the CMB; D pp is the particle diffusion coefficient due to the resonant scattering of Alfven 
waves that are assumed to be isotopically distributed [23]: 
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Numerical solutions for the set of coupled differential equations given in Equations 2.17 and 
2.18 have been presented in e.g., [30] and [38]. Here, however, aiming to study the cumulative 
effect of turbulent acceleration of cluster mergers, we simplify the calculation of individual 
events by using an analytical steady-state solution following [23] and [24], 


In the steady-state regime, we further simplify Equations 2.17 and 2.18 by setting the con¬ 
dition dWk/dt = 0 and assuming that the wave and particle spectra are simple power-laws: 

W fc = W fc ° k~ v (2.20) 

f(p) = foP~ s (2.21) 


Noting that the lower limits of Equations 2.16 and 2.19 are given by p m ; n 
kmin = Om e /p respectively, the power-law assumptions lead to: 
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As demonstrated in [23, 24], if and only if 


v = 3 


f hn e /k and 


( 2 . 22 ) 
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and 

a self-similar solution with df /dt 


S = sA, 

0 can be reached : 
f{jp) oc p ~ (6 ~ St) 


(2.26) 

(2.27) 


With this spectrum and using Equation 2.11, the energy loss rate of electrons is found to 
satisfy the relation 


^si 


J Sp 2 cf (p)Anp 2 dp 


= Pa 


^(6 — s t )(4-s t ) 


(2.28) 


Taking m = 5/3 for Kolmogorov turbulence and using Equations 2.9, 2.27, 2.28, as well 
as the pa factor below equation 2.10, we find a synchrotron and inverse Compton emission 
power that is comparable to the energy flux injected in the Alfven waves Psi = Pa , and an 

electron energy spectrum N(E ) oc EN 2 - 5 . Specifically, Equation 2.28 normalizes the electron 

1 /2 

flux through the relation Pa ~ (167r/3)5/oPmax, and thus the equilibrium spectrum can be 
written as: 


_dNe_\ _ 3 Pa c 5 /2 

dE e dV) inj 4S(P max )V2 e 


(2.29) 


The spectrum of synchrotron emission is given by: 

dN s = V3e 3 B r E ““ dE E s dN e 

dE s dtdV hm e c 2 E s J me e hv c dE e dV 


(2.30) 
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Figure 1. Properties of the turbulence generated during the merger of two halos as a function 
of the halo mass ratio. Top panels: Turbulent velocity at injection (left) and the injection scale of 
the turbulent Eddies (right). Bottom panels: Timescale of the turbulence ta = lo/v t (left), Ratio 
between the turbulent energy and the thermal energy of the host halo (right). Values are evaluated 
at two redshifts: z = 0.1 (dashed lines) and z = 1 (solid lines), and for three different masses of the 
host halo: Ml = 10 15 M© (red), 10 14 M 0 (blue), and 10 13 M 0 (green). r] v = 0.8 has been assumed 
for all calculations. 

where F(x) = x J a °° K 5 / 3 (x')dx' is the Synchrotron function, K 5 / 3 is the modified Bessel 
function, and v c is the critical frequency of synchrotron emission, v c = 3y 2 eB /( 47 rm e c) = 
1.6 (B/lpG)(E e /10GeV) 2 GHz. Note that the synchrotron emission function F{x) peaks at 
x = v/v c ~ 0.29. 

While .E/max in Equation 2.15 represents the maximum energy an electron can gain through 
damping, our self-similar solution is strictly time-independent. This indicates that the 
acceleration timescale for electrons t a cc ~ p 2 /2D PP is equal to the energy loss timescale 
floss ~ E/{dE/dt) = c/(SE) at all energies. Thus, from the analytic formula derived here 
it is difficult to directly compute the maximum electron energy after energy losses are taken 
into account. However, the relevant timescales have been computed in [30], who find that 
it is possible to accelerate electrons to at least 7 ~ 10 16 , and possibly much higher as the 
calculations of [30] scan a limited range of magnetic field strengths and time-ranges. While 
additional magnetohydrodynamic studies are necessary to fully explore these models in the 
context of the ARCADE excess, we note that our results are not highly sensitive to the max¬ 
imum electron energy. In Figure 7 of the appendix, we will show that the fit to the ARCADE 
data is only mildly impacted at high-frequencies when the maximum electron energy is de- 
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creased. In part, this is due to the fact that a significant fraction of the total synchrotron 
signal is produced by clusters with intermediate masses, as we will show in Figure 5), which 
have a corresponding 7 max much smaller than 10 5 . Finally, we note that Alfven reaccelera¬ 
tion can also produce high-energy electrons with energy 7 > 10 5 through a second process. 
Protons accelerate to very high energy via Alfven reacceleration can produce high-energy 
electrons through inelastic collisions with the interstellar medium [38]. This process poten¬ 
tially provides a new source of high-energy electrons capable of increasing the synchrotron 
emission above GHz. 

We stress that calculations in the above two sections are entirely based on the Eilek-Henriksen 
solution [23], which is strictly true only under ideal circumstances, when waves mainly in¬ 
teract with electrons (but not protons), wave-electron interaction dominates over wave-wave 
cascading, and energy losses follow p 2 scaling (like inverse-Compton and synchrotron pro¬ 
cesses). A magnetohydrodynamic approach like in [30, 38] would be needed to fully describe 
the system. Yet the scenarios discussed in [30, 38] could be further complicated by the fact 
that Alfven waves can develop a scale-dependent anisotropy during their cascading (e.g., 
[39-41]), which could strongly reduce the acceleration rate [42]; while on the other hand, 
compressive turbulence in the ICM could enhance the acceleration (for example fast modes, 
see [43, 44]). The current work, however, aims to explore the possible link between the 
ARCADE-2 excess and reacceleration of electrons induced by cluster mergers. Therefore we 
limit our calculation of the electron spectrum to the ideal case. 


2.3 Density Profile, Magnetic Field and Temperature of a Cluster 


We assume a simple /3-model for the density profile of the thermal gas in the ICM [38, 45]: 


Picm(^) = Picm,o 



(2.31) 


where j3 = 0.8 and r c is the core radius which is taken as 10% of the cluster’s virial radius 
r v i r [24]. The virial radius is determined through the relation M = 47 rr^ ir A v i r (z) p m (0)/3, 
with A v j r (’:) £3 (187T 2 + 82x — 39x 2 )/(x + 1) and x = Qm{z) — 1 [46]. The density profile is 
normalized by: 

/icmM = J piCM(r)4:TTr 2 dr (2.32) 

with /icm ~ 0.13(M/10 14 M q ) 0 - 16 the baryon fraction of clusters [47]. 


For clusters, which we define as any halo with a total mass M > 10 12 Mq, we adopt a 
magnetic field distribution as a function of the radial distance from the center of the cluster 


to be [48-50] 


B(M,r) = B 0 




(2.33) 


where ij ~ 3/8 [14], and Bq is the central magnetic field strength of clusters. We leave Bq 
as our first adjustable parameter with a default value Bo = 10/xG. The central magnetic 
field strengths should also depend on the temperature, and therefore the mass of clusters 
as suggested in e.g. [51, 52], Here we assume that the field strength scales to the mass of 
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cluster by (M/M*)“, with M* = 10 14 M 0 and a 0.3 (inferred from the estimations of held 
strengths of Coma and A3667 in [51]). For galaxy-sized halos with magnetic fields dominated 
by their dense baryonic cores, we assume a magnetic held model [13] 

B(r) = max (.Bi e~ r/R \B 2 e" r/i?2 ) (2.34) 

where Bi = 7.6 fiG and Ri = 0.025 r v ; r and E >2 = 35 fiG with R 2 = 0.008r v ; r . 


As we will show in Section 3, the cumulative radio signal is dominated by clusters with 
mass greater than 1O 14 M 0 . Therefore our result is mostly independent on the magnetic 
held models for galaxy-sized halos. As for the magnetic held model for clusters, our default 
assumption of Bq = 10 /iG is close to the upper limit of cluster magnetic held strength 
derived from Faraday rotation measurements [17]. Nevertheless, a significant enhancement 
of the magnetic held would be expected due to the boost in the gas density by the shocks 
during merger events [53, 54], For comparison, our default model predicts B v i r = 2.5/uG 
at the virial radius of a 10 15 M 0 halo during the merger, while observations of radio relics 
suggest (i G-level magnetic held at the position of the shock fronts [55] . 


The temperature of ICM gas can be described by [24] 

Ficm = F vir + -T g 


(2.35) 


where 2>T g /2 = 0.8 keV is the temperature of the preshock gas [24] and T v i r is the virial 
temperature of a cluster, 

_ GM [ApTHp . 

vir “ 2 r vil k B { 

with Up = 0.61 being the mean molecular weight and Jzb being the Boltzman constant. 


2.4 Cluster Mergers 


Turbulence is expected to be injected during the passage of subhalos [56]. In what follows, 
we denote the smaller (sub-) halo with a subscript ’S’, and the larger (host) halo with the 
subscript ’L’. The relative velocity of the two halos with masses of Ml and Ms can be 
described by [57] 


Vi = 


2G{M l + M s ) 
Rl 



(2.37) 


with rj d 

by 


4 ((Ml + Ms)/M l) 1 / 3 ■ Following [24] and [56], we estimate the turbulent velocity 

Vt = Vv Vi (2.38) 


with 0 < r) v < 1 being our second free parameter. The top left panel of Figure 1 shows vt 
as a function of the halo mass ratio Ml/Ms, taking ij v = 0.8. The turbulent velocity at 
the injection scale is around 300 — 2000kms _1 depending on the mass of the host halo, a 
value that is consistent with previous studies (e.g. [58]). Note that when Ml/Ms 1, the 
turbulent velocity is approximately the free-fall velocity due to the gravitational potential of 
M l , yielding v t ~ GM L /R L . 
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Figure 2. Energy power of the injected turbulence (dashed), the Alfven waves generated by the 
turbulence (dash-dotted), and the synchrotron emission from electrons accelerated by the Alfven waves 
(solid). As in Figure 1, values are calculated for three host halo masses, Ml = 10 13 ,10 14 ,10 15 M 0 
in green, blue and red, respectively. Parameters assumed in this plot include central magnetic field 
strength B 0 = 10 fiG and rj v = 0.8. 


The smaller, infalling halo suffers ram-pressure stripping when it crosses the larger one. At 
the stripping radius, the ram pressure applied to the smaller halo is comparable to its static 
pressure: 


PLvf = 


Ps(r s )k B T s 
fir rip 


(2.39) 


Turbulence is expected to be injected on a similar scale as the stripping radius. In [24] and 
[30], lo = Vi r s was assumed, taking r]i to be another free parameter. However, noting that r/i 
is poorly constraint observationally, and that the Alfven power is significantly less dependent 
on Iq than on vt (recalling that Pa oc l Q 7 ^), in this work we reduce the degrees of 
freedom in the model and set Iq = 2 r s . We show Iq as a function of the halo mass ratio in 
the top right panel of Figure 1. Note that in the extreme case Ml Mg, r s —>• 0, meaning 
that the smaller halo is almost immediately stripped when it approaches the larger cluster. 


The duration for the injection of turbulence is of the order of the crossing time, ta = lo/vt 
[24]. As plotted in the bottom left panel of Figure 1, the typical timescale for turbulence 
ranges from 0.1 to a few Gyrs depending on the ratio of the halo masses. Note that ta is 
longer for smaller host halos, due to their significantly slower vt- The volume swept by the 
passage of the subhalos would be filled with turbulence. High-resolution hydrodynamical 
simulations found that in major mergers, mergers dominate turbulent diffusion in the entire 
virial volume; in minor mergers, several effects triggered by the merger including sloshing, 
rotational large-scale motions and hydrodynamical instabilities would introduce turbulence 
to a size close to the scale radius of the subcluster [59, 60]. We therefore assume that during 
ta, Alfven waves can be excited within a swept volume Va ~ ttR^Iq. Then in the bottom 
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Figure 3. The merger rate for a host halo of mass Ml with a smaller halo Ms, as a function of the 
halo mass ratio M^/Mg in the comoving frame. Legends are the same as in Figure 1. 

right panel of Figure 1 we plot the ratio between the turbulent energy Et = Pl v 2 Va and 
the thermal energy of the two halos E t h = {Ml + Ms) ks TjJ fim p . For 10 14 — 10 15 M & 
halos this ratio is found to be 10 — 18%, consistent with numerical simulations [61] as well 
as observations of the Coma cluster, which suggested a lower limit of 10% [62]. 

As a consistency check, Figure 2 compares the power contained in the turbulence injected by 
the merger with the Alfven waves generated by the turbulence, and the synchrotron emission 
from electrons accelerated by the Alfven waves. Note that the power is evaluated over the 
turbulent injection time, and is thus slightly higher for a larger halo mass ratio which has a 
shorter ta- 

The probability that the larger cluster Ml undergoes a merger with a smaller cluster with 
mass Mg per unit time is given by [57] 


d 2 p 

[2 5 c {z) 

dln5 c {z ) 

dlna(M') 

dMgdt V 

7r a{M') 

dt 

dM 1 


(2.40) 
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where M' = Ml + Mg, 8 c (z) is the critical overdensity given by the expression: 
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Figure 4. The contribution to the isotropic radio background from electrons accelerated by cluster 
turbulence. The synchrotron emission is shown to be comparable to the temperature excess reported 
by the ARCADE-2 collaboration [5] . Results with three sets of parameters are presented in blue, red 
and green, corresponding to central magnetic field strength Bq = 10,4 pG, and a ratio between the 
turbulent velocity and the initial relative velocity of the merging clusters being p v = 0.86, 0.74. 


for a flat universe, + flv = 1, and a (M, z) is the mis variance of the linear density field, 

a 2 = J dkP hn (k)W(kR)k 2 (2.42) 

where W(kR) is the Fourier transform of the real-space top-hat window function of radius 
R = (SM/Airpm) 1 / 3 [63] with p m (z) = p m ( 0) (1 + z) 3 , and P lin is the linear matter power 
spectrum [64]. The merger rates of three host halo masses Ml = 10 13 ,10,10 15 Mq at 
z = 0.1,1 are shown in Figure 3. 


2.5 Cumulative Synchrotron Intensity 


The mean intensity of synchrotron emission from the turbulent acceleration in merging halos 
can be written as 


I(Es) 



W[M l ,(1 + z)E s ,z] 


(2.43) 


where x is the commoving distance, dn/dM is the halo mass distribution function which is 
given by 


dn _ . p m din a 1 

dM ~ dM 

and for f(cr) we use the mass function multiplicity described by [65]. 


(2.44) 
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Figure 5. The sub-contributions to the total ARCADE-2 temperature stemming from host halos 
of various masses. The magnetic field is assumed to have a central value B 0 = 10 /iG. The turbulent 
velocity is assumed to be a fraction rj v = 0.86 of the initial relative velocity of the merging clusters. 
In this figure we show the contributions from three mass groups: M < 10 12 Mq, 10 12 Mq < M < 
1O 14 M 0 , and M > 10 14 Mq. The M > 10 14 Mq group (in red) dominates the signal (lying almost 
directly on top of the solid blue line). 


The window function W (Ml, E s ) describes the synchrotron flux expected from all the merger 
events of a halo Ml at redshift z, 


W (Ml, E s , z) = 


x 


Hz r d 2 p 

— - / dMs -rr~z —— 

47r A/s, min dMsdt 


dN s 

dE s dt 


t a (M l , M s , z) 


(2.45) 


where dN s /dE s dtTA sums the synchrotron emission from all radii of a larger halo during its 
merger with a smaller halo, with 


dN s 

dE s dt 



dr Airr 2 


dN s 

dE s dVdt 


(B,n pi T, l 0 ,v t ) 


(2.46) 


In the calculation we set the lower limit of Ml to be ML, m in = 10 11 Mq. This lower limit 
is however irrelevant, as we will show in Sec 3 that the diffusive emission is dominated by 
massive clusters with M > 1O 14 M0. The lower limit of Mg is rather important, because 
for the same halo, its merger rate with smaller halos is significantly higher than with larger 
halos, as shown in Figure 3. Considering that most of the signal comes from clusters with 
M > 10 14 Mq, and that the stripping radius of a 10 6 Mq halo merging into a 10 14 Mq halo 
is already close to 0, we set Mg )m ; n = 10 6 M 0 . 
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Figure 6. The sub-contributions to the total ARCADE-2 temperature stemming from halos at 
various distances. As in Figure 5, B 0 = 10 fiG and rj v = 0.86 have been assumed for calculation. In 
this figure we show the contributions from four redshift bins: z < 0.1, 0.1 <z<l, l<z<2, and 
2<z<5. 

3 Results 

In Figure 4, we present the extragalactic radio background produced from electrons acceler¬ 
ated by turbulence in galaxy clusters as a function of frequency. As described in Section 2, 
our model has two adjustable parameters: (1) the magnetic field strength at the center of 
galaxy clusters (defined as halos with mass above 10 12 M 0 ), Bo, and (2) the ratio between the 
turbulent velocity and the initial relative velocity of the merging clusters, r] v . Results with 
two sets of parameters are shown in Figure 4, including Bo = 10,4/xG and r] v = 0.86,0.74. 
In particular, the values of r] v have been chosen to fit the data point at 0.05 GHz. We find 
that (Bo = 10 fiG, r] v = 0.86) gives a better fit to the ARCADE-2 excess, while the other 
case with weaker central magnetic field provides a even softer spectra. 

The softer synchrotron spectrum in models with lower Bo is due to the constraint that the 
magnetic held strength at the virial radius of galaxy clusters is proportional to the central 
value (as given by Equation 2.33). Weaker magnetic held channel a larger fraction of the 
turbulent energy into synchrotron emission at lower frequencies. For comparison, in a cluster 
with M = 10 14 M 0 at z = 0.1, the magnetic held strength at the virial radius would be 
B v ir =1.2/rG for Bq = 10 /j>G, corresponding to a synchrotron peak frequency u p = 2 GHz for 
a 10 GeV electron, but would be B V1V =0.5 fiG and u p =0.8 GHz for Bq = 4 /iG. 

In Figure 5 we decompose the cumulative radio emission into sub-contributions from host 
halos in three mass groups: M < 10 12 M 0 , 10 12 M 0 < M < 10 14 M 0 , and M > 1O 14 M 0 . 
We find that the most massive group contributes dominantly with 60%, the 10 12 M 0 < 
M < 10 14 M & group contributes 40%, and the least massive group contributes negligibly to 
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the observed radio excess. One caveat is that the fractional contribution from clusters in 
different mass groups relies on the scaling of magnetic field with the cluster mass (as defined 
in Equation 2.33) through the the R-dependent E max (Equation 2.15 and 2.26). The partition 
of contributions of subgroups in Figure 5 could be altered if very different E m ax and B(M) 
relations would be invoked. 

Additionally, in Figure 6 we decompose the synchrotron signals into sub-contributions from 
host halos located at four different redshift bins: z < 0.1, 0.1 < z < 1, 1 < z < 2, and 
2 < z < 5. It is found that the merger events at 0.1 < z < 1 contribute most significantly, as 
a result of the fact that massive clusters are most abundant in the nearby universe. 

Intriguingly, we note that the anisotropy produced by cluster mergers is also expected to be 
quite low. While we withhold on a quantitative calculation of the anisotropy, we note that a 
10 14 Mq cluster at z = 1 has an angular size of 6', which corresponds to i ~ 3600. Therefore 
the dominate emission sources in our model (Mx > 10 14 Mq and z < 1) have angular sizes 
significantly greater than the 2' constraint imposed by [16], and are beyond the range of the 
angular scales that could be probed by the radio measurements quoted in [15]. 

4 Discussion and Conclusions 

In this paper we have shown that cluster mergers are capable of powering the bright isotropic 
synchrotron emission necessary to explain the ARCADE-2 excess. Specifically, these cluster 
mergers provide a number of advantages over current explanations from both astrophysical 
and dark matter models. The large power of the ARCADE-2 excess is easily provided by 
cluster mergers, and the spectrum of the ARCADE-2 excess is well-fit for steady-state electron 
spectra following dN/dE oc E ~ 2 ' 5 with a cut off at E max ~ 50 GeV and a cluster magnetic 
field strength near the virial radius of ~1 fiG. It is enticing to note that the observed radio 
data is so easily fit by these parameter choices, which are among the more theoretically well- 
motivated and observationally constrained cluster parameters [54, 55]. The turbulence and 
subsequent re-acceleration of non-thermal electrons is intrinsically a non-thermal process, 
which produces negligible IR and X-Ray emission, thus remaining consistent with high- 
energy constraints. Finally, the synchrotron emission is dominated by the largest clusters, 
significantly decreasing the predicted anisotropy from this model at high-multipoles. 

There are two caveats to the analysis shown in this work. First, while our model predicts 
a maximum electron energy of ~50 GeV, this stems from a simple, time-independent ap¬ 
proximation that does not carefully take into account between Alfven reacceleration and 
energy-loss process such as synchrotron emission. A full magnetohydrodynamic model of 
cluster mergers is necessary to carefully model the electron energy spectrum in both minor 
and major mergers, an effort which lies beyond the scope of the current work. In Appendix A, 
we discuss this caveat in more detail, and show that reasonable variations in the maximum 
electron energy moderately impact our fit to the high-frequency excess, but do not effect 
fits below ~1 GHz. Additionally, we have assumed that the power of Alfven waves could 
be efficiently converted into the re-acceleration of electrons. This assumption may not be 
guaranteed. For instance, when the energy density of relativistic protons is larger than that 
of electrons, an efficient electron acceleration would unavoidably accelerate protons, leading 
to an increase of the damping rate of the Alfven waves and thus a decrease in the efficiency 


16 


of electron acceleration. However, as cosmic rays at GeV-TeV energies are expected to be 
confined within the cluster volume for cosmological times, re-energization of secondary elec¬ 
trons and positrons is possible due to the inelastic interactions between cosmic ray protons 
and the intracluster medium [38]. 

Furthermore, in our model we set the turbulence scale to be 2 times the stripping radius, 
which could underestimate Iq especially in cases of minor mergers. For example, observations 
show evidence of minor merging activity in the massive cluster Abell 2142 (M 200 = 1.3 x 
10 15 Mq) with a subhalo of mass ~ 10 13 M 0 [66, 67]. This merging activity is suggested to 
be inducing turbulence and creating an extended radio halo with a linear scale of 2 Mpc 
[68], but in our model such a merger only generates turbulence for less than 200 kpc. The 
possibility of larger turbulent volume leaves room for a scenario with less efficient electron 
acceleration. On the contrary, if the actual turbulence volume is not as large as modelled 
in this work, or equivalently, if on average the turbulence channels less than 20% of thermal 
energy in major mergers, and 0.5% in minor mergers, our model would have over estimated 
the synchrotron emission from cluster mergers. 

Finally, we note that our model provides a highly testable observational signature. Mea¬ 
surements of the radio anisotropy at larger angular scales would be capable of detecting a 
significant anisotropy stemming from cluster contributions to the ARCADE-2 excess. Forth¬ 
coming observations with the Low Frequency Array (LOFAR) [69], the Murchison Widefield 
Array (MWA) [70], and the planned Square Kilometre Array (SKA) [71] should be sufficient 
to detect such anisotropy signals [72, 73]. Moreover, the cross-correlation of this anisotropy 
measurement with existing catalogs of massive clusters (e.g. SPT catalog [74], ACT survey 
[75, 76], Planck survey [77], ROSAT all-sky survey [78]), would provide smoking gun evidence 
for cluster contributions to the isotropic radio background. Due to the expected presence of 
high degrees of turbulence in cluster shocks, the lack of significant radio emission in cluster 
mergers would also be of great interest, informing our understanding of re-acceleration in 
these energetic environments, and motivating the search for more exotic explanations of the 
ARCADE-2 data. 

A Impact of more realistic electron spectra 

Throughout this work we have adopted an analytical approach to calculate the Alfven re¬ 
acceleration following [23] and [24], In particular, we have assumed that injection electrons 
follow a simple power-law of dN/dE e oc F% 2 - 5 (see equation 2.29). However, numerical 
simulations suggested a more comprehensive electron spectrum with a characteristic bump 
slightly before the cut-off energy (see for example, Figure 3 of [38]), as a result of the interplay 
of the turbulent re-acceleration of electrons and protons, as well as the inelastic interaction 
between cosmic ray protons and the ICM. To evaluate the impact of more realistic electron 
spectra, below we test our model by adding artificial bumps to the injection spectrum to 
mimic the numerical results. Specifically, we assume that a bump with dN/dE e N times 
higher than the original power-law flux is placed between E max /ci and E max /c r . The modified 
spectrum is then normalized by the Alfven injection power through equation 2.28. Figure 7 
present the synchrotron emission from the following six tests: 1) N = 5, q = 10 3 , c r = 10 2 
2) N = 100, q = 10 3 , c r = 10 2 3) N = 5, q = 10 2 , c r = 10 4) N = 100, q = 10 2 , 
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Figure 7. Comparison of radio emission by different electron spectra. In addition to the power- 
law spectrum with dN/dE e oc E~ 2 5 in our default case (see equation 2.29), we assume the pres¬ 
ence of a bump with an amplitude N times higher for the following energy ranges: 1) N = 5 at 
[£ ; max/10 3 ,£max/10 2 ], 2) N = 100 at [E max /10 3 , A ma x/10 2 ] , 3) N = 5 at |£ max /10 2 ,.E max /10] and 
4) N = 100 at [A max /10 2 , -Emax/10]. We also show the impact to the fit if the maximum electron 
energy is 5 times greater or smaller than the A max defined in equation 2.15. The modified spectra 
are normalized by the Alfven injection power through equation 2.28. In all cases the magnetic field 
is assumed to have a central value B 0 = 10 fiG. The turbulent velocity is assumed to be a fraction 
7j v = 0.86 of the initial relative velocity of the merging clusters. 


c r = 10. The spectrum shows a significant derivation from the initial result in case 2), where 
most electrons have low energies and thus produce a much softer spectrum compared to 
observation. However for all the other cases, after integrated over the entire halo population, 
the bumps are smoothed out and produce very small changes to our results. The above tests 
show that our conclusions are robust under different inputs of electron spectra. 

In addition, to test the dependence of our results on the maximum energy of the re-accelerated 
electrons, we artificially injected spectra with maximum energies 5 times higher or lower 
than the E max calculated by equation 2.15. As shown by the black and yellow lines in 
Figure 7, E max could slightly modify the overall synchrotron spectrum, but the effect is 
rather minor. On the other hand, if E max would have been significantly less than the value 
from equation 2.15, the resulting spectrum would be even softer, and it may become difficult 
to explain the high-frequency tail of the ARCADE-2 emission. A more advanced MHD 
modelling would be important to test the impact from a time-dependent E mSLX . 


B Magnetic energy compared to the thermal energy 


As a consistency check, in Figure 8 we estimate the ratio of the magnetic energy, Eb to the 
total thermal energy E t h of the merging clusters in the turbulent volume. Specifically, the 
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Figure 8. The ratio of the magnetic energy Eb to the total thermal energy E t h in the turbulent 
volume as a function of the mass ratios, for a central magnetic field strength B 0 = 4 fj,G (solid 
lines), and B 0 = 10/iG (dashed lines). For reference, the black dash-dotted line indicates the level 
Eth/Es ~ 40 as suggested in [79]. 


ratio is calculated by 


Eb 

Eth 


( M l ,M s ,z ) 


B 2 JTr \ ,f(M L + M s )k B T L V t \ 


( 0 . 1 ) 


The ratios reported in figure 8 range from 2% to 8% for Bo = 10 fiG, and from 0.2% to 1.1% 
for Bo = 4/iG. These estimated values are comparable with the level (~ 2.5%) suggested in 
[79]. 
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